Search for: All records

Repeat RNA sequences self-associate to form condensates. Simulations of a coarse-grained single-interaction site model for (CAG)n (n = 30 and 31) show that the salt-dependent free energy gap, ΔGS, between the ground (perfect hairpin) and the excited state (slipped hairpin (SH) with one CAG overhang) of the monomer for (n even) is the primary factor that determines the rates and yield of self-assembly. For odd n, the free energy (GS) of the ground state, which is an SH, is used to predict the self-association kinetics. As the monovalent salt concentration, CS, increases, ΔGS and GS increase, which decreases the rates of dimer formation. In contrast, ΔGS for shuffled sequences, with the same length and sequence composition as (CAG)31, is larger, which suppresses their propensities to aggregate. Although demonstrated explicitly for (CAG) polymers, the finding of inverse correlation between the free energy gap and RNA aggregation is general. 

Abstract Folding of ribozymes into well-defined tertiary structures usually requires divalent cations. How Mg2+ ions direct the folding kinetics has been a long-standing unsolved problem because experiments cannot detect the positions and dynamics of ions. To address this problem, we used molecular simulations to dissect the folding kinetics of the Azoarcus ribozyme by monitoring the path each molecule takes to reach the folded state. We quantitatively establish that Mg2+ binding to specific sites, coupled with counter-ion release of monovalent cations, stimulate the formation of secondary and tertiary structures, leading to diverse pathways that include direct rapid folding and trapping in misfolded structures. In some molecules, key tertiary structural elements form when Mg2+ ions bind to specific RNA sites at the earliest stages of the folding, leading to specific collapse and rapid folding. In others, the formation of non-native base pairs, whose rearrangement is needed to reach the folded state, is the rate-limiting step. Escape from energetic traps, driven by thermal fluctuations, occurs readily. In contrast, the transition to the native state from long-lived topologically trapped native-like metastable states is extremely slow. Specific collapse and formation of energetically or topologically frustrated states occur early in the assembly process. 

Low-complexity nucleotide repeat sequences, which are implicated in several neurological disorders, undergo liquid–liquid phase separation (LLPS) provided the number of repeat units, n , exceeds a critical value. Here, we establish a link between the folding landscapes of the monomers of trinucleotide repeats and their propensity to self-associate. Simulations using a coarse-grained Self-Organized Polymer (SOP) model for (CAG) n repeats in monovalent salt solutions reproduce experimentally measured melting temperatures, which are available only for small n . By extending the simulations to large n , we show that the free-energy gap, ΔG S , between the ground state (GS) and slipped hairpin (SH) states is a predictor of aggregation propensity. The GS for even n is a perfect hairpin (PH), whereas it is a SH when n is odd. The value of ΔG S (zero for odd n ) is larger for even n than for odd n . As a result, the rate of dimer formation is slower in (CAG) 30 relative to (CAG) 31 , thus linking ΔG S to RNA–RNA association. The yield of the dimer decreases dramatically, compared to the wild type, in mutant sequences in which the population of the SH decreases substantially. Association between RNA chains is preceded by a transition to the SH even if the GS is a PH. The finding that the excitation spectrum—which depends on the exact sequence, n , and ionic conditions—is a predictor of self-association should also hold for other RNAs (mRNA for example) that undergo LLPS. 

Both the small and large subunits of the ribosome, the molecular machine that synthesizes proteins, are complexes of ribosomal RNAs (rRNAs) and a number of proteins. In bacteria, the small subunit has a single 16S rRNA whose folding is the first step in its assembly. The central domain of the 16S rRNA folds independently, driven either by Mg²⁺ions or by interaction with ribosomal proteins. To provide a quantitative description of ion-induced folding of the ∼350-nucleotide rRNA, we carried out extensive coarse-grained molecular simulations spanning Mg²⁺concentration between 0 and 30 mM. The Mg²⁺dependence of the radius of gyration shows that globally the rRNA folds cooperatively. Surprisingly, various structural elements order at different Mg²⁺concentrations, indicative of the heterogeneous assembly even within a single domain of the rRNA. Binding of Mg²⁺ions is highly specific, with successive ion condensation resulting in nucleation of tertiary structures. We also predict the Mg²⁺-dependent protection factors, measurable in hydroxyl radical footprinting experiments, which corroborate the specificity of Mg²⁺-induced folding. The simulations, which agree quantitatively with several experiments on the folding of a three-way junction, show that its folding is preceded by formation of other tertiary contacts in the central junction. Our work provides a starting point in simulating the early events in the assembly of the small subunit of the ribosome. 

RNA molecules cannot fold in the absence of counterions. Experiments are typically performed in the presence of monovalent and divalent cations. How to treat the impact of a solution containing a mixture of both ion types on RNA folding has remained a challenging problem for decades. By exploiting the large concentration difference between divalent and monovalent ions used in experiments, we develop a theory based on the reference interaction site model (RISM), which allows us to treat divalent cations explicitly while keeping the implicit screening effect due to monovalent ions. Our theory captures both the inner shell and outer shell coordination of divalent cations to phosphate groups, which we demonstrate is crucial for an accurate calculation of RNA folding thermodynamics. The RISM theory for ion–phosphate interactions when combined with simulations based on a transferable coarse-grained model allows us to predict accurately the folding of several RNA molecules in a mixture containing monovalent and divalent ions. The calculated folding free energies and ion-preferential coefficients for RNA molecules (pseudoknots, a fragment of the rRNA, and the aptamer domain of the adenine riboswitch) are in excellent agreement with experiments over a wide range of monovalent and divalent ion concentrations. Because the theory is general, it can be readily used to investigate ion and sequence effects on DNA properties.